Dependence of resistivity on surface profile in 
nanoscale metal films and wires 



Baruch Feldman, Rui Deng and Scott T. Dunham 
Physics Department, University of Washington, Seattle, WA 
Electrical Engineering Department, University of Washington, Seattle, WA 

We extend quantum models of nanowire surface scattering to incorporate bulk resistivity and 
extract an expression for the increased resistivity due to surface roughness. To learn how to im- 
prove conductivity, we calculate conductivity degradation from individual wavelengths of surface 
roughness, and show how these can be convolved to give resistivity for arbitrary surfaces. We re- 
view measurements from Cu films and conclude that roughness at short wavelengths (less than 100 
nm) dominates scattering, and that primarily specular scattering should be achievable for RMS 
roughness below about 0.7 nm. 



As the minimum feature size in semiconductor technol- 
ogy continues to shrink, metal nanowires with thickness 
<45 nm are now needed to interconnect electronic nan- 
odevices. However, measurements show nanowires have 
substantially higher resistivity than bulk metals [l|, [3], 
leading to interconnect delays, power loss, and other lim- 
its on performance. Scattering from surfaces, interfaces 
and grain boundaries are the causes of this conductiv- 
ity degradation, but microscopic understanding of these 
effects and quantitative predictions of their magnitude 
have been limited. Here, we investigate the detailed de- 
pendence of conductivity on surface roughness profile and 
analyze the resulting technological impact. 

The first quantitative treatments of surface and size ef- 
fects in thin films or wires were the semiclassical methods 
of Fuchs Q and Sondheimer These approaches as- 
sume a ratio p of carrier collisions with the surface reflect 
specularly, while X—p scatter diffusely. Such theories can 
be fit to experiment with p as a free parameter, but do 
not provide insight into how to improve conductivity. 

More recently, surface roughness scattering has raised 
the attention of researchers in industry [3, [3, 0]> an d 
quantum mechanical approaches to surface scattering 
calculations have been proposed. The two primary ap- 
proaches include the Kubo linear response theory of 
Tesanovic et al. [3] and Trivedi and Ashcroft [3], and 
the diagrammatic Keldysh formalism of Meyerovich and 
collaborators 13, II , G3] • Here we follow the approach 
of Meyerovich et al, which is readily applied to arbitrary 
surface roughness profiles. We calculate the contribution 
of each spatial frequency of surface roughness and con- 
volve with roughness data extracted from experiments to 
gain insight into the nature of surface roughness scatter- 
ing. 

In our conductivity calculations, we consider a thin 
film because it reproduces the major qualitative results 
of a wire (and matches quantitatively when Eq. @ be- 
low holds), while avoiding strong localization and other 
effects that make ID systems problematic to deal with 
theoretically 13, 13] • For the technologically important 
10-100 nm scale, wire conductivity can be accurately es- 



timated by combining effects of scattering from sidewalls 
to that from top and bottom surfaces. 

In a thin film of thickness L, boundary conditions at 
the surfaces lead to a density of states quantized in the 
transverse direction. As a result, the conduction band, 
described as the set of states at the Fermi energy, is bro- 
ken into subbands with continuous parallel and quantized 
transverse components of the Bloch wavevector. Conduc- 
tion states are then described by a subband index j and 
a 2D wavevector kj, subject to the constraint that the 
total energy is equal to the Fermi energy: 
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(We treat the Fermi surface as effectively spherical, which 
is particularly appropriate for the best conducting met- 
als, Ag, Cu, and Au). Even in a perfectly smooth film, 
this quantization leads to thickness-dependent conduc- 
tivity, and to the quantum size effect (QSE), caused 
by the quantized dependence of the density of states 
on thickness [3] which is significant for very thin films 
(<5nm). 

Theoretical approaches to rough surfaces [3, E3] em- 
ploy a (non-unitary) transformation to map the film 
with position-dependent surface into a flat film with bulk 
(non-Hermitian) perturbations. The scattering depends 
on the power spectral density (PSD) of the roughness, 
defined as the Fourier transformed surface height corre- 
lation function. In [13, [13] , a general isotropic 2D rough- 
ness power spectrum C(l^l) is treated with diagrammatic 
perturbation theory. By Fermi's Golden Rule, the spatial 
frequencies of roughness determine the interband transi- 
tion rates and hence a momentum loss rate matrix: 
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Here Tikj is the in-plane momentum satisfying (TTJ) for 
subband J , x is the angle between initial and final carrier 
momentum, ( (kj,kji ,x) = C(y^f + — 2kjky cos x), 
and superscripts denote (2D) angular harmonics: 
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Surface roughness-limited conductivity is given by 
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where p s is resistivity, the scalar t s is the overall sur- 
face relaxation or mean free time, n = k F /3ir 2 is carrier 
density, and l/r s oc p s is the overall surface momentum 
loss rate. Note that our definition differs from [l2| by 
an extra factor of 3/2irL because we use the usual 3D 
conductivity, as in [9(. 

To combine bulk and surface scattering, we extend 
the method in [9(, adding momentum loss rates within 
subbands, to the case with interband transitions by 
adding matrices: r _1 = t^ 1 + T.T 1 . Since the primary 
bulk scattering mechanism at room temperature, acous- 
tic phonons, is nearly isotropic [20j . we use 

with vp = 1.6 x 10 8 cm/s the Fermi velocity and A& = 
39 nm the bulk mean free path for copper. 

Adding matrices produces very different results from 
adding the scalars 1/r oc p. Matthiessen's rule, which 
states that pr = Pi + Pi for independent (series) resis- 
tivity mechanisms, breaks down in thin films when com- 
bining bulk and surface scattering 0, 0, [ll[ (Fig. [1] be- 
low). This breakdown can be understood because, ab- 
sent bulk scattering, conductivity is dominated by carri- 
ers with momentum nearly parallel to the surface (low j) 
which rarely scatter from the surface. To consider surface 
scattering together with bulk scattering, we define 
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the effective surface roughness contribution to resistivity, 
which is independent of pb to first order [25j |. 

To study the effect of individual spatial frequencies of 
roughness on resistivity, we perform a first order func- 
tional expansion on p e s s . We define the first variation in 
p c s s with respect to the PSD component at wavevector 
with magnitude ko as the response to a special PSD: 
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Here p^? ^ is the response to a 2D PSD of the form 

_ l 2 5(\k\ - fc ) 
Cfeo(k) = k ' 



with 2tt1 2 the mean squared roughness of this PSD. The 
factor (l 2 /ko) in ([7]) is necessary for consistent units. 

Consistent with the validity of ([2]) - (f5]) to first order in 
roughness, we perform a first order functional expansion 
of p c s s in (: 
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Here the LHS is the resistivity from an arbitrary 2D 
isotropic PSD (, and p°^ k ^ in the RHS is given by Eq. 0. 
The angular harmonics for (|8|) are given by 
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where 9 is the Heaviside step function, and the delta 
function sets the angle x between the initial and final 
wavevectors 
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Our results for the functional derivative ([7]) for surface- 
only scattering l/r s oc p s and effective rate with bulk 
scattering l/r° oc p° are plotted in Figure[T] We tested 
the relation ([9]) for several PSDs, including our fit to 
experimental surface roughness data (see below). For 
the experimental fit, we find (J9j) holds to 2% for a 25 
nm film and 1% for 100 nm. We conclude that, for films 
thick enough that p e s <C p, Eq. provides an effective 
calculation of resistivities. 

For larger k values, which dominate surface scatter- 
ing, the thickness dependence of the surface only resistiv- 
ity is p s oc L~ 2 , as in [^, [l2j]. But for large L, we expect 
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FIG. 1: The momentum loss rate l/r s oc p s and l/rf 1 oc p'f 1 
in response to a single-frequency PSD (0, (j8]) in copper with 
I = 10 pm as a function of fcrj. Note drastic breakdown of 
(8) Matthiessen's rule, with l/r 3 cft > l/r s . 
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p to approach the semiclassical result of Fuchs [3j and 
Sondhcimer Wl: 
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Fig. [T] indeed shows p^ s oc k^/L for most of the fco range. 

For very small fco values, pf^ oc ko/L 2 , while the 
surface-only rate actually increases with L, p s oc L. This 
can be explained by quantization. The subbands are 
spaced closest for lowest j, so for fco < \/37r/ (HfL), there 
is no scattering between higher order subbands. As fco be- 
comes smaller, interband scattering due to surface rough- 
ness becomes impossible for lower and lower subbands. 
For p s , an increase in L decreases the spacing between 
kj values, introducing interband transitions between low 
j states where the spacing is smallest. Physically, this 
couples states traveling nearly parallel to the film sur- 
face with other states that interact much more strongly 
with the surface, providing a mechanism to increase p s 
with L. For p c s R , in contrast, electrons in low j states 
are already frequently scattered by bulk scattering. For 
small fco, only intraband scattering is possible, so there 
are always two final states and p^f oc ko/L 2 . At higher 
fco, the number of available final subbands becomes pro- 
portional to L, so p° s oc k 2 /L. 

As can be seen from Fig. [U spatial frequencies near 
(shortest wavelengths) have the strongest momentum 
loss, but the impact on conductivity depends on the ac- 
tual roughness PSD of metal films and wires. Any at- 
tempts to improve conductivity will benefit from a knowl- 
edge of which components of surface roughness give the 
most improvement for the resources spent. Theoreti- 
cal calculations often assume a Gaussian roughness spec- 
trum, but experiments show that many PSD forms are 
present depending on the wire deposition conditions [13] , 
and that real PSDs can fall of more slowly than Gaus- 
sian [l8j |. Unfortunately, the experimental literature on 
surface roughness spectra for metals is limited and fo- 
cuses on relatively large length scales. Thus, we look to 
other materials. Feenstra et al. [l8| observed that ID 
STM scans of InAs/GaSb superlattice interfaces showed 
Lorentzian distributions, 
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For isotropic roughness, this corresponds to a 2D PSD of 
the form 
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with mean squared roughness 27rA 2 . Eq. (fTTjl also fits 
the AFM results of Moseler et al. for Cu films with 
A = 18 nm, A = 1.8 A, as shown in Fig. [2 Other 
experiments on copper films confirm a correlation length 
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FIG. 2: Normalized roughness spectra for Cu films [19( com- 
pared to fits to Eq. with correlation length A = 18 nm 
and Gaussian with A = 1/a — 25 nm. Both forms fit the 
available data, but give very different predictions of fall-off in 
roughness at higher spatial frequencies. 



of ~20 nm [21]. We also fit [19('s data to a Gaussian 
PSD, as shown in the figure. Because measurements at 
high spatial frequencies (which have a particularly strong 
effect on scattering) are lacking, our goal is to extrapolate 
from these fits [261 ] - 

We can use extrapolations from the Moseler data to 
calculate the resistivity, or equivalently the specular frac- 
tion p. Substituting I© in ([IT 
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which (for L >100 nm) is independent of L. We find p 
values of essentially 1 for the Gaussian PSD and p — 0.96 
(1 — p = 0.04) for Eq. (fTT|) . We get the same results when 
we use the full PSDs directly as in Eqs. © - © 

The fact that our analysis predicts highly specular 
(p ~ 1) surface scattering for technologically- achievable 
surface roughness suggests that surface roughness scat- 
tering is a surmountable barrier to high conductivity. We 
can extract the most important components of roughness, 
taking into account both the relative strength of scatter- 
ing and the observed roughness spectra. The effective 
diffuse scattering rate as a function of spatial frequency 
is shown in Fig. [3] More accurate measurements of the 
high frequency portion of the PSD are clearly needed, 
as the frequencies above 0.1 nm -1 are most critical to 
conductivity degradation. 

Another way to understand these results is to note that 
for (fTTj) with A = 18 nm, p = 90% corresponds to an 
RMS roughness of 7 A, compared to exp erimental mea- 



surements in the range 2 - 11 A 

The experimental literature is somewhat mixed on the 
relative importance of surface scattering. Many results 
suggest that the observed resistivity increase is domi- 
nated by grain boundary rather than surface scatter- 
ing [2j], while some extract values of p near (diffuse 
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FIG. 3: Integrand of Eq. |12J) with PSDs from Fig.H Diffuse 
fraction (1 — p) oc p c s s is the area under the (ID) curve as 
shown. The tail of the PSD is seen to have a major effect. 
The most important roughness components are 0.1 nm -1 
k < k F = 13.6 nm -1 for Eq. dTTJ and ko ~ 0.05 nm' 1 
Gaussian. 



< 
for 



rather than specular scattering) [231 ] . As we have seen, 
p depends strongly (0(l 2 )) on RMS roughness, which in 
turn depends on anneal times, deposition conditions, and 
other process variables. Another interesting explanation 
for this discrepancy may come from the experiment of 



Rossnagel et al. [24(, who found that conductivity de- 



creased strongly upon the deposition of an ultrathin Ta 
layer on top of a Cu film, but that the conductivity re- 
covered when the Ta film was exposed to air, thereby 
oxidizing to become insulating. These observations sug- 
gest that thin barrier/ adhesion layers rather than sur- 
face/interface roughness may be causing the apparent 
diffuse surface scattering. 

To summarize, we have found that bulk scattering can 
be included in quantum models of surface scattering by 
adding t _1 matrices. This leads to violation of Math- 
iessen's rule, but an effective surface resistivity p° inde- 
pendent of bulk scattering can be extracted. The resis- 
tivity from individual wavelengths of roughness can be 
convolved with roughness PSD to get p e s s for arbitrary 
surface. Our analysis suggests that roughness with wave- 
length within 1-2 orders of magnitude of the Fermi wave- 
length is the most critical for conductivity degradation. 
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We can show by expanding in ^-t^ 1 that I/t^ = 
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the highest frequency is due to the hnite sampling inter- 
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